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Abstract — Decoding sparse quantum codes can be accom- 
plished by syndrome-based decoding using a belief propagation 
(BP) algorithm. We significantly improve this decoding scheme by 
developing a new feedback adjustment strategy for the standard 
BP algorithm. In our feedback procedure, we exploit much of 
the information from stabilizers, not just the syndrome but also 
the values of the frustrated checks on individual qubits of the 
code and the channel model. Furthermore we show that our 
decoding algorithm is superior to belief propagation algorithms 
using only the syndrome in the feedback procedure for all cases 
of the depolarizing channel. Our algorithm does not increase the 
measurement overhead compared to the previous method, as the 
extra information comes for free from the requisite stabilizer 
measurements. 

Index Terms — Sparse quantum codes, quantum error correc- 
tion, quantum channels, belief propagation, stabilizers. 



I. Introduction 

COMMUNICATION is limited by noise in channels, 
but error correction methods can efficiently offset this 
restriction in both classical Q], and quantum O, flU, 
0, |6), H cases. At a simple level, multiple copies of the 
information can be transmitted, and a majority rule can be 
applied to discern the correct code, but such coding is neither 
practical nor efficient. Sparse graph coding, such as Gallager's 
low-density parity-check (LDPC) codes, offers an efficient 
alternative that approaches the Shannon information limit (T), 
12. Fortunately quantum coding and decoding strategies can 
be constructed from their classical counterparts, but unfortu- 
nately this mapping from classical to quantum coding can 
be problematic due to the requirement that quantum codes 
satisfy the duality-containing condition 0, 10. Moreover, 
due to increased challenges posed by these quantum codes, 
performance improvement requires further progress in the 
proposed decoding algorithm IflOl . IfTTII . 

The entanglement-assisted (EA) stabilizer formalism adds 
error-free entangled bits as a consumable resource for perform- 
ing quantum error correction. This EA approach overcomes 
the duality-containing requirement and thus offers a rich lode 
of quantum error correction protocols inspired by classical 
protocols fl2l . Using the EA approach, modern codes for 
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classical channels, such as sparse codes, can easily be imported 
as quantum error-correcting codes lfl3l . fl4l . 

Our aim is to improve belief propagation (BP) decoding 
methods so that quantum coding is dramatically improved over 
existing techniques. Specifically our numerical results pre- 
sented here show that an improved BP method whose heuris- 
tical feedback strategies based on exploiting all accessible 
information from stabilizer measurements, yield a dramatically 
improved block error rate (BER) for any depolarizing channel. 
Our methods should work for any Pauli noise channel. 

For n qubits, a Pauli channel is defined by the mapping 

S(C 2 ") 3pn^ p E EpE\ p E > 0, J] p E = 1 (1) 

EeC„ 

with G„ the n-fold tensor product of single-qubit Pauli opera- 
tors S e {I,X,Z, Y = XZ\. Our interest is focused on memory- 
less channels wherein the error on each qubit is independent 
of the error on any other qubit. In particular we consider the 
depolarizing channel, which is the most-studied case J9), IflOl . 
ifTTI . ifTJl . l27l : for fixed channel error probability p, the error 
on qubit q is given by 

S (C 2 ) 3 Pq h» PlPq + (1 - pi) (x Pq X + Y Pq Y + Z Pq Z) /3 (2) 

with pi : = 1 —p denoting the probability of no error occurring. 

In quantum settings, due to the inability to measure each 
and every qubit, syndrome-based decoding is typically chosen. 
Consequently, decoders using Pauli channels are generally 
considered to be hard-decision decoders. In other words, in 
quantum decoding, the conventional soft-decision techniques 
are not applicable when the channel is in the Pauli channel 
model. Generally, sparse quantum codes are decoded by using 
syndrome-based BP decoding algorithms which automatically 
imply hard-decision decoding. Though there is an equivalence 
between syndrome-based decoding and a posteriori probabil- 
ity decoding (signal-based decoding) under this setting, the 
syndrome-based decoding results in a serious drawback to BP 
decoder: the symmetric degeneracy error IflOl . Fortunately, 
Poulin and Chung (PC08) [TOj propose a solution to the 
symmetric degeneracy error by using the random perturbation 
method. 

Motivated by the soft-decision techniques used in classical 
settings, which provide extra reliable information on the 
message nodes thereby yielding a better error correcting ca- 
pacity lfl6l . we develop a new heuristical feedback adjustment 
strategy for the standard BP decoder. When used in decoding 
sparse quantum codes, our method can on one hand solve the 
symmetric degeneracy error problem, and on the other hand 
can provide more useful information to the message nodes. 
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The difference between PC08 and our approach is that, 
in PC08, they feed back only the syndrome of the decoder 
output to adjust prior error probability distributions for re- 
ceived qubits. These adjusted distributions are then fed back 
into the decoder. We significantly improve this protocol by 
feeding back not just the syndrome but also the values of 
the frustrated checks on individual qubits of the code and the 
channel model; accordingly we introduce a new adjustment 
strategy. Specifically, our approach, which is based not only 
on syndromes but also on frustrated checks obtained from full 
stabilizer measurements and the channel model, yields a better 
BER for the case of depolarizing quantum channels. 

We provide a detailed description of our basic BP decoder 
for decoding sparse quantum codes. The basic BP decoder 
introduced here is inspired by the strategy used for decoding 
sparse classical quaternary codes under the BP algorithm. 
Using this strategy, we can decode sparse quantum codes 
directly regardless of whether they arise from classical binary 
codes or not and regardless of whether they are Calderbank- 
Shor-Steane (CSS) construction Q, |[19) codes or not. 

II. BP Iterative Decoding 

In this section, we briefly reprise the essential elements 
of BP iterative decoding for sparse quantum codes. In Sub- 
sec. IH-AI we discuss the key idea of standard BP iterative 
algorithms. Then we compare decoding of classical codes 
vs quantum codes, and introduce the standard BP decoding 
for quantum codes. Finally in Subsec. IH-BI we show how to 
decode sparse quantum codes in GF(4) based on the standard 
BP algorithm. 

A. BP decoding algorithm for sparse quantum codes 

Consider a k-bii message that is encoded into an n-bit 
codeword, which is then transmitted through a noisy channel. 
The received message is an error-prone vector over the output 
alphabet, and there is no guarantee that decoding will reveal 
the original codeword. However, the codeword can be guessed 
with a high probability of being correct by maximizing the 
probability for this codeword based on the observed output 
vector 11201 . Unfortunately for a linear block code, encoding k 
information bits into n bits allows 2 k possible n-bit codewords, 
and calculating conditional probabilities for individual code- 
words is too expensive in practice. 

BP algorithms overcome this inefficiency for sparse 
codes [T|. The strategy is to represent a linear-block classical 
error-correcting code (CECC) by a Tanner graph comprising 
message nodes and check nodes corresponding respectively to 
received bits and check constraints. Then an iterative algorithm 
recovers the values of received bits as follows 1201 . ETI . At 
each round, probabilities are passed from message nodes {v} to 
check nodes |c) and then from |c) back to {v}. The probabilities 
from {v} to {c} are computed based on the observed value 
of the message node and other probabilities passed from the 
neighboring check nodes to the respective message node. It 
turns out that, for decoding sparse codes, the BP algorithm 
can provide a reasonable trade-off between complexity and 
performance. 



For a linear block CECC, the code space can be viewed as 
the orthogonal projection space (solution space) of its check 
matrix. The sender Alice transmits her message as a codeword 
that is encoded according to a specific check (or generator) 
matrix through the channel. When receiver Bob obtains the 
channel output, the received vector may not be the solution 
vector of the check matrix. Therefore Bob needs to apply 
a smart algorithm to recover the codeword efficiently. As 
the check matrix is sparse, Bob employs the BP algorithm. 
First, Bob measures each bit to obtain its posterior probability 
distribution. Subsequently, he puts these probabilities into the 
BP decoder, and, based on the constraint that the codeword 
should be the orthogonal vector of the check matrix, he infers 
the original message. This procedure extends naturally to 
sparse quantum codes. 

In the quantum case, the stabilizer formalism for a quantum 
error-correcting code (QECC) is useful. The code space within 
C 2 " corresponds to the simultaneous +1 eigenspace of all m 
generators of an abelian subgroup J5) 

S =(Si SJcG». (3) 

Alice transmits her message as a codeword, which propagates 
through the depolarizing quantum channel. Due to channel 
errors, the received codeword may not be the simultaneous 
+ 1 eigenspace of S . Bob measures S to obtain the syndrom^ 
s = {s\,...,s m ) e {-1,1} and then estimates which error E 
from Eq. (HJ has occurred based on the noise model and 
the syndrome, which is effectively a check constraint. Bob 
redresses the error by applying the same error operation he 
inferred onto his received codeword to try to recover the 
quantum code. For Pauli channels all errors are n-fold tensor 
products of single-qubit Pauli operators hence square to the 
identity; thus re-applying the same error should restore the 
original codeword. 

Due to the close analogy between CECCs and QECCs as 
discussed above, it is tempting to extend the BP iterative 
decoding algorithm to quantum settings. In fact a QECC can 
also be represented by a decorated Tanner graph lfl7l . Ifl8l . 
Therefore, it is natural to transmit the messages from {c} to {v} 
and the messages from the {vj to {c} according to the principles 
of BP algorithms ifTTl. |[T8l. 

Here we first briefly review the BP algorithm used for 
decoding sparse quantum codes iflOl but begin with necessary 
terminology. The probability of Pauli error E q occurring on 
qubit q is p q (E q ). The neighborhood of qubit q is denoted 
by n(q), and we define n(c) similarly as the set of qubits 
connected to check S c . The tensor product of Pauli errors 
over n(c) is denoted E c . Then, following PC08's notation, for 
two 1-qubit case, E-E' = 1 means [E,E r ] = and E-E' = -1 
means \E,E'} = 0. For two n-qubit Pauli operators such as 
E — Ei ® Ei ® • • • <8> E„ and E' = E[ ® E' 2 ® • • • ® E' n , we have 

// 

E E' :=Y] E k -E' k . (4) 

k=\ 

'if [S C ,E] = 0, s c = 1 or else .s c = — 1, where c = l,-- ,m and S c can 
be viewed as a check node in the Tanner graph of a QECC with stabilizer S 
from Eq. {3). 
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The messages from check nodes to qubit nodes are denoted 

m c -> q (E q ) oc ^ (s Sc ,s c -e c mf^ciEq,)), (5) 

Eqt q'en{c)\q 
q'en(c)\q 

which is defined only up to a constant factor. The factor can be 
fixed by normalization [lOj. The messages from qubit nodes 
to check nodes are similarly denoted as 

m q ^ c {E q ) oc p q (E q ) m c ,^ q {E q ), (6) 

c'en(q)\c 

where p q (E q ) is the initial probability in the definition of the 
memoryless channel. Then the beliefs b q (E q ) are constructed 
by first initializing m q ^ c (E q ) = p q (E q ), then evaluating ac- 
cording to 

b q (E q ) = Pq (E q ) I m c ^ q {E q ), (7) 

cen(q) 

after the iteration procedure based on Eqs. (0 and ([6j. In fact, 
Eqs. (0 and (0 define a sum-product iterative procedure for 
decoding sparse quantum codes. Hence this algorithm is also 
called the sum-product algorithm (SPA), which is one of the 
most important algorithms based on BP. 

In order to clarify the BP decoding algorithm for sparse 
quantum codes, we now show how to implement Eqs. (0 and 
(0 in GF(4). 

B. Decoding sparse quantum codes in GF(4) 

There is a convenient isomorphism between the Pauli group 
d generated by {I,X,Z, Y = XZ\ and the Galois Field GF(4) 
generated by {0, l,a>, a) = a> 2 }. The isomorphism is explained 
by the element identification 

/ ^ 0, X ^ 1, ,Z ^ to, Y <-> to(to 2 ) (8) 

and the operation identification multiplication^addition and 
commutativity<->trace inner product [0. For the one-qubit 
case, we have, 

P,QeGi <-> P,£>eGF(4), 

[P, Q] = <=> Tr(P X Q) = 0, [P, Q] = <f> Tr(P X £>) = 1 

(9) 

The isomorphism is readily extended to the n-qubit case: 

P,QeG„ <-> u P ,v Q e GF(4)", 

[P, Q] = «* Tr(H P • v Q ) = 0, {P, Q) = «* Tr(H P • v ) = 1. 

(10) 

where '■' used here (between two vectors) is a regular 
inner product. That is, for u = (ui,U2,---,u n ) and v = 
(vi , V2, ■ • • , v„), we have, 

k=n 

h-v:=^«|X vi. (11) 

k=\ 

The addition and multiplication rules of GF(4) are shown in 
Tables J] and [TTJ respectively. 

We can import the strategy used for decoding sparse clas- 
sical quaternary codes under BP to use for decoding sparse 
quantum codes. This adaptation to quantum codes is achieved 
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by transforming the check criterion (0 from commutativity to 
trace inner product. 

For example, suppose s c = 1, Then the check criterion for 
this check node should be S c ■ E c = 1, which is equivalent to 
[S C ,E C ] = 0. According to isomorphism ( fTOb , Tr(M £i -v^) = 0, 
which implies that either u E ■ v Sc = or u E ■ v Sc = 1 because 
Tr(0) = and Tr(l) = 0. We define p q [x] as the probability 
for E q x5 c to take the value x, with S c the q th entry of S C Q 
and p[x] the probability of u Ec • Vs c = (excluding the q th 
entry) to take the value x. As q e n(c), 

S Cq e {X, Z, Y) ^> § Cq € {\,<jj,<b},x € {0, \,a>,a>}, (12) 

Thus, 

P,[0] = (p[0]+p[l])/2 = p ? [l], (13) 

Pq [co] = (p[co] + p[u>]) 12 = p q m. (14) 

Similarly, if s c — — 1, then {Sc^d = 0, which implies 
that either u Ec • Vs c = to or u E ■ Vs c = to because Tr(a>) = 1 
and Tr(w) = 1 . According to Table [I] 

p q [0] = (p[o)]+p[cb])/2 = p q [ll (15) 

p q [co]=(p[0]+pm)/2 = p q m (16) 

With this knowledge, the message m c _> 9 can be computed by: 

m c ^q(E q ) = x 5 C? ]. (17) 

Given m c ^ q (E q ), m q ^ c (E q ) can then be derived directly by 
substituting Eq. (flTt into Eq. (0. 

Decoding based on GF(4) affords the advantage that sparse 
quantum codes can be decoded without losing correlations be- 
tween errors that would otherwise impact the coding scheme's 
performance. In particular sparse quantum codes can be de- 
coded directly regardless of whether they arise from a classical 
binary codes or not and regardless of whether they follow a 
CSS construction code or not. For comparison, we henceforth 
refer to the decoder based on this BP decoding algorithm as 
the 'standard BP decoder'. 

2 For E q and S c the mapped elements of E q and S c in GF(4), see Eq. (9). 
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III. Enhanced feedback BP iterative decoding for sparse 

QUANTUM CODES 

In 2008 PC08 proposed a random perturbation strategy 
based on the syndrome of the BP decoder output to overcome 
the symmetric degeneracy problem. The random perturbation 
method is simple and efficient but unfortunately overlooks 
many errors. As pointed out in PC08, all errors in their 
simulations could be attributed to the decoder rather than to 
the finite minimal distance of the code. 

We observe the same problems as PC08 do when using 
the random perturbation method. In particular, most errors in 
the simulations are detected errors, not the undetected errors. 
This result suggests that feeding back useful information to 
the standard BP decode could help. The feedback would help 
the BP decoder determine a valid output E ut, whose syndrome 
s(E out ) is identical with the observed syndrome s. Then E out 
will equal E with high probability, for E the error occurring 
on the transmitted quantum state during transmission. In 
this section, we present our enhanced feedback BP iterative 
decoding algorithm, which provides useful information to the 
BP decoder based on exploiting not only the syndrome but 
also the stabilizer itself and the channel model in the feedback 
procedure. 

Then suppose that, subsequent to the standard BP decoding 
procedure, s(£out) is not identical with the observed syn- 
drome s. Then evidently at least one entry of E out occurring 
on one of the qubits connected to check S c has an error, for 
S c a frustrated check. If we can reset the initial probability 
distributions for errors occurring on the qubits connected to 
S c in a more reliable way, but not simply by adopting the prior 
probability distributions obtained from the channel model, 
more useful information could be provided to the BP decoder. 
Furthermore symmetric degeneracy could also be overcome as 
well. Thus it is reasonable to expect that BP decoding ability 
will be significantly improved. 

Our new feedback BP iterative decoder is devised as fol- 
lows. First the standard BP decoder is used to identify the er- 
ror E based on the observed syndrome and the prior probability 
distributions of each entry of E. When the standard BP decoder 
fails (s(E out ) + s), we randomly choose a frustrated check, 
say S c . There are only two results for this frustrated check: 
s c = — 1 and S c ■ E out = 1 or else s c = 1 and S c ■ E out = - 1 . 
Either way we choose a random qubit that connects to S c , for 
example qubit q, and then reset p q (E q ) according to the value 
of S Cq (the entry to which qubit q corresponds in S c ) and the 
channel model as follows. 

Let Si 2,3 £ {X,Z, Y = XZ) be three distinct elements. Sup- 
pose 5 Ci| = Ei, and the channel is the depolarizing channel 
with crossover probability p = 1 — pj, for pj established in 
Eq. ©. Then, if s = -1 and S c ■ E out = 1, we reset 

Pq (I) = (1 - Pl )/2 = p ? (E 2 ) = Pl /2 = p ? (E 3 ). (18) 

3 In addition to the random perturbation technique, PC08 also proposed 
freezing and collision techniques. However, their simulations and our own 
simulations indicate that the random perturbation technique has the best 
performance of the three approaches. Therefore, we adopt the random pertur- 
bation technique as the comparator. 



For example, if 5 Cl| = X, we obtain 

p q (l) = (1 - Pi)/2 = Pq (X), Pq (Z) = Pl /2 = Pq (Y). (19) 
Instead suppose s c = 1 and S c ■ E out = — 1. Then we reset 

p q (I) = Pl /2 = p 9 (S0, Pq (Z 2 ) = (1 - P1 )I2 = p ? (E 3 ). (20) 
In this case, if S Q = X, we obtain 

Pq (I) = Pl /2 = Pq (X), Pq (Z) = (1 - P1 )I2 = p q (Y). (21) 

Next, following a similar strategy to PC08's, we feed back 
this adjustment to the standard BP decoder and let it iterate 
with this modified probability distribution for T pelt (a pre- 
determined fixed number of the reiteration) steps. If BP halts 
during this period (i.e. s(E 0Ut ) corresponds to the syndrome s), 
the procedure is complete and terminates. 

If c is still frustrated, we restore the probability distribu- 
tion p q (E q ) and choose a different qubit connected to c, say 
qubit q', and reset p' q (E' q ) based on the same strategy that has 
already been used for resetting p q (E q ). If S c is not frustrated, 
but the halting condition is not yet satisfied, we choose another 
frustrated check S c > and adjust the probability distribution of 
the error occurring on one of its qubits as we do for a frustrated 
check S c . 

Now suppose that the decoder still cannot yield a valid 
output after trying a pre-determined number n a feedback 
adjustments for qubits according to the procedures mentioned 
above. In other words, suppose that, after traversing n a entries 
associated with frustrated checks, no valid result is yielded by 
the decoder. In this case, we allow the decoder to terminate 
this procedure and report a failure. Without the termination 
condition, the decoder is wasting time traversing all entries 
associated with frustrated checks and thereby introducing 
unwanted redundancy into the procedure. 

Eqs. (fT8l and (1201 are justified for the following reason. 
Empirically, the decoder is naturally too biased towards the 
identity / fTOl due to the prior probability distribution obtained 
from the channel model. This biased probability distribution 
will lead to S c ■ E out = 1 while s c = -1. Heuristically, we 
reset p q (E q ) to make it anti-commute with 5 Ci| with a greater 
probability, since at least one entry of E anti-commutes with 
its corresponding entry in S c - The channel is depolarizing that 
is Pauli errors X, Z, Y occur with equal probability, so, in 
Eq. (IT8b . we let E2 and E3 share the probability as /□ 

When s c = 1, setting E q to / is conducive to obtain- 
ing S c ■ E out = 1. Another Pauli operator could also be 
conducive to making S c ■ E out — 1 but is ignored by the prior 
probability distribution. We can readily ascertain this Pauli 
operator according to the entries of S c . According to Eq. i20\ . 
we let this Pauli operator (Si) share the same probability pj 
with / and have £2 and S3 share the small probability 1 - pj 
equally. 

Our strategy retains PC08's capability of breaking the sym- 
metric degeneracy while feeding back more useful information 
to the BP decoder than PC08 in order to help the decoder 
determine a valid output. More useful information comes from 
retaining full information about the entire stabilizer rather 

4 Because (1-p/) <K pi, we let / and Zi share (1-p/) equally for simplicity. 
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than just a syndrome. Specifically, our feedback adjustment 
strategy is not solely based on the syndrome but also on the 
channel model and on the individual values of the entries of 
the frustrated checks (the checks are just the generators of the 
stabilizer). 

IV. A Case Study 

In this section we illustrate our enhanced feedback iterative 
decoding algorithm to show how it succeeds in helping the 
BP decoder find a valid output. For [[«,£; c]] denoting an 
EA QECC that encodes k qubits into n qubits with the 
help of c ancillary ebits Ifl2ll . we construct a simple EA 
QECC for [[4, 1; 1]] as follows. According to the EA QECC 
formalism, we start with a classical [4, 2] quaternary code with 
check matrix 

1 co 1 
110 1 



(22) 



First we transform H c to 



H' = 



(jjH c 



( 1 OJ 
1 1 

a) a> 



1 
1 
d 



<x> CO 



CO 



(23) 



then transform H' Q to a set of (perhaps non-commuting) gen- 
erators 

X Z X I \ 
XXIX 
Z Y Z I 
Z Z I Z 



(24) 




Number of iterations, n 
(C) 



Number of iterations, n 
(D) 




4 8 12 16 

Number of iterations, n 



12 16 20 

Number of iterations, n 



Fig. 1 . Beliefs as a function of number of iterations n for an EA QECC with 
stabilizer generators XZXIX, XXIXZ, YZZXI and ZXXYI and with syndrome 
(-1, + 1, +1, +1) using the standard BP decoding algorithm. The prior p q (E q ) 
of the four sent qubits are obtained from the depolarizing channel with p = 
0.1. (A), (B), (C) and (D) show the beliefs for qubits 1,2,3 and 4, respectively. 
Qubit 5 is assumed to be held by the receiver and thus to be error-free. The 
symbols □, A, O and V indicate the beliefs for /, X, Z and Y, respectively. 
In this case, the output of the standard BP decoder is E out = IYII, which is 
a detected error. 



Next we transform H q into a canonical form H q by mul- 
tiplying the third generator by the second generator and 
by multiplying the fourth generator by the first and second 
generator. Thus, we obtain 



( X Z X I \ 
XXIX 
Y Z Z X 
Z X X Y 



(25) 



Now, it is easy to check that the first generator anti-commutes 
with the second generator, and the last two generators com- 
mute with each other as well as commuting with the first 
two generators. According to the EA QECC formalism, this 
coding scheme needs one ebit to assist encoding. The extended 
commuting set of generators is {H q } for 



X 
X 
Y 
Z 



z 

X 

z 

X 



'1' 

X I 

I X 

Z X 

X Y 



X 

z 
I 
I 



(26) 



which is just the stabilizer of [[4, 1; 1]]. 
It is easy to check that H q is isomorphic to 



X 

z 
I 
I 



I 
I 

z 
I 



I I 

I I 

I I 

Z I 



X 

z 
I 
I 



(27) 



which means we can encode one logical qubit into four 
physical bits with the help of one entangled pair plus two 



ancilla bits via a unitary transformation 11121 . As the half of 
the entangled pair being held by the receiver is assumed to 
be error-free, the error-correcting capacity of [[4, 1; 1]] only 
depends on the left part of H q , which corresponds to the four 
transmitted qubits held by the sender. 

Suppose the qubits of [[4, 1; 1]] are sent through a depo- 
larizing channel and error E - IIZX is applied to the four 
transmitted qubits. Measurement of the stabilizer reveals the 
error syndrome s = (-1,+1,+1,+1). In Fig. Q] we display the 
performance of the standard BP decoder under this setting. 
From Fig.Q] it is easy to check that the syndrome of the output 
of the standard BP decoder is (-1, -1, -1, -1)0 We now use 
the standard BP decoder with PC08's random perturbation. 
Comparing between the error syndrome s = (-1,+1,+1,+1) 
and the syndrome of the output of the standard BP decoder 
(-1,-1,-1,-1), the probability distributions of the errors 
occurring on the qubits connected to the second, third or fourth 
checks can be reset. For all these cases, the decoder could not 
yield an appropriate recovery. Here as an example, we use 
PC08's random perturbation strategy for the frustrated check 
S 2 and show the corresponding performance in Fig. [2] 

Finally, we use our enhanced feedback BP decoding algo- 
rithm. By the same method, our feedback strategy can reset the 
probability distributions of the errors occurring on the qubits 
connected to the second, third or the fourth checks. When the 
decoder chooses the fourth entry (X) of the frustrated check S 2 

5 It is not surprising that the standard BP decoder has failed because 
[[4, 1; 1]] is designed to correct one error but not two errors. 
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Fig. 2. Same as for Fig. [T] but using the standard BP decoding algorithm 
replaced by PC08 and applying random perturbation strength S = 1 to the 
prior for qubits 1 , 2 and 4. 
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Number of iterations, n Number of iterations, n 

(C) (D) 




3 4 8 12 16 2( 3 4 8 12 16 20 

Number of iterations, n Number of iterations, n 



Fig. 3. Same as for Fig.[2]but using our enhanced feedback iterative decoding 
algorithm. Here we reset P4(E$) according to Eq. )20t . In this case, only three 
iterations were required to yield a valid output £ oul = IIZX, which is exactly 
the error occurring on [[4, 1;1]] during its transmission. 

and resets P4CE4) according to Eq. i20i . the correct decoding 
result arises in just a few iterations. We show the performance 
of our approach in Fig. [3] 

In this example, PC08's method could not help the standard 
BP decoder to yield a valid output, However, our enhanced 
feedback BP decoding algorithm yields a valid output in just 
a few iterations. In fact even if the error occurring on the 
first four qubits is not IIZX but YZII, which has the same 



syndrome s — (-1,+1,+1,+1) as IIZX, our decoding output 
IIZX can still recover the transmitted quantum state because 

IIZXI x YZIII = YZZXI, (28) 

which is just the third generator of the extended stabilizer. 
Because it is hard to check whether E ■ E out e S when the 
number of generators of S is large, we choose E = E out as the 
success criterion of the decoding result in our simulations. 

V. Simulation Results 

We have applied our enhanced feedback iterative BP de- 
coding algorithm to a variety of sparse quantum codes, in- 
cluding conventional sparse quantum codes and EA sparse 
quantum codes over depolarizing channels. In each case, our 
improved BP decoder yields significantly lower BER over 
both the standard BP decoder and the BP decoder with 
PC08's random perturbation. In the following subsections, we 
simulate decoding of sparse quantum codes with different code 
parameters (including block lengths, rates, and row weight) 
under the three decoders to demonstrate the superiority of our 
approach. 



A. Example: Conventional Sparse Quantum Codes 

Conventional sparse quantum codes can be constructed 
from sparse classical dual-containing codes. One of the most 
successful dual-containing constructions is the so-called "Con- 
struction B" J5J, which is built as follows. First we take an 
n/2 x n/2 cyclic matrix C with row weight L/2, and define 

H Q = [C,C T ]- (29) 

Then we delete some rows from Hq to obtain a matrix 
H with m rows. By construction, H is dual-containing. A 
conventional sparse quantum code with length n can thus be 
given according to CSS construction. 

Our first example is based on this strategy. We first construct 
a cyclic binary LDPC code [63, 37] with row weight 8 based 
on finite geometries jl26l . This code has a 63 x 63 cyclic 
check matrix with 26 independent rows and 37 redundant 
rows. Hence, we can construct a conventional sparse quantum 
code with n = 126, m = 26, and L = 16. We refer to this 
[[126,74]] quantum LDPC code with rate k/n = 74/126 and 
row weight 16 as 'Conventional"]^ We show the performance 
of the standard BP decoder, the BP decoder with PC08's 
random perturbation and our enhanced feedback BP decoder 
when applied to Conventional over depolarizing channels in 
Fig.gO 

6 Note that in order to improve the performance of this code, we have 
preserved those redundant rows. 

7 In our simulation, we take the maximum number of entries traversed by 
the enhanced feedback decoder round down to one fifth of the code length for 
short codes, one tenth of the code length for medium codes and one fortieth 
for long codes. 
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A— Standard BP decoder 
^PC08 
v — Enhanced Feedback 



10" 



10 

0.015 0.02 0.03 0.04 0.05 

Cross-over probability of depolarizing channel, p 

Fig. 4. Block error rate vs cross-over probability of depolarizing channel p 
for three different BP decoders when applied to Conventional. The maximum 
number of iterations is 90. The maximal number of iterations between each 
perturbation is 40. The strength of the random perturbation used here is 6 = 
0.1. The maximal number of entries traversed by enhanced feedback decoder 
round down to one fifth of the code length, which equals 25 here. 



B. Examples: EA Sparse Quantum Codes 

EA sparse quantum codes can be constructed from arbitrary 
binary or quaternary sparse classical codes 1121 . Here we 
present the decoding of three EA sparse quantum codes 
that are constructed in the same way as the construction 
of [[4, 1; 1]] discussed in Sec. [IV] These three codes have dif- 
ferent net rate (including positive, zero and negative net rates), 
different block lengths and different row weight distributions. 

The first EA sparse quantum code we consider is 
[[128, 79; 49]], which is constructed from an irregular qua- 
ternary sparse classical code [128,79] with girth g > 6. The 
classical code [128,79] is designed by using the progressive- 
edge-growth (PEG) algorithm under the constraint that the 
degree distribution of bit nodes is |[T4l 



(2 : 33%, 3 : 33%, 4 : 22%, 5 : 22%} 



(30) 



As there are 49 pairs of anti-commuting generators in H q 
according to the construction method of Eq. (|25T >. then the 
EA formalism yields the EA sparse quantum code 

[[n,2k-n + c;c]] = [[128, 2 x 79 - 128 + 49; 49]], (31) 

which can be constructed from the classical quaternary code 
[n,k] = [128,79] HI. The code [[128, 79; 49]] has a positive 
net rate 

^ = Z^i? (32) 
n 128 

Here we refer to this code as "EA-1" and show the per- 
formance of the three decoders with respect to EA-1 over 
depolarizing channels. The performance, in terms of block 
error rate, is summarized in Fig. [5] 

Our second EA sparse quantum code is [[8 16, 404; 404]], 
which is constructed from a binary regular sparse classical 
code [816,408]. The classical [816,408], which has a rate 0.5, 
row weight 6 and column weight 3, was constructed by 
Mac Kay Ir22ll . As there are c = 404 pairs of anti-commuting 




Standard BP decoder 
PC08 

Enhanced Feedback 



0.02 0.03 0.04 0.05 0.06 

Cross-over probability of depolarizing channel, p 

Fig. 5. Same as for Fig. 4 but replacing Conventional by EA-1. 



generators and 8 commuting generators in H q , we obtain an 
EA sparse quantum code [[8 16, 404; 404]] according to the 
EA stabilizer formalism. Here we take this code as our second 
example of EA sparse quantum codes whose net rate is 

k- c 404 - 404 



816 



0, 



(33) 



and call it "EA-2". We show the performance of the three 
decoders when applied to EA-2 over depolarizing channels in 
Fig.|6]in terms of BErU 




Standard BP decoder 
-^^PC08 

-V — Enhanced Feedback 



0.09 0.10 0.11 0.12 0.13 

Cross-over probability of depolarizing channel, p 

Fig. 6. Same as for Fig. 4 but replacing Conventional by EA-2 and changing 
the maximum number of entries traversed by the enhanced feedback decoder 
round down to one tenth of the code length, which equals 81 here. 

The third EA sparse quantum code is [[1920,638; 1278]], 
which is constructed from [1920,640], a famous sparse clas- 
sical code proposed as the international standard for cellular 
telephones. The [1920,640] code, which has a rate 1/3 and 
row weight 4, was also constructed by Mac Kay 1231 . Fol- 
lowing the construction procedure mentioned in Sec. IIV1 the 
number of anti-commuting pairs in H q is 1278; thus, we obtain 

8 Net rates for EA quantum codes can be negative or zero. See Appendix 
for more detail. 
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an EA sparse quantum code [[1920,638; 1278]]. Here we take 
this code as our third example of EA sparse quantum codes, 
this time with a negative net rate 

k-c 638 - 1278 1 

= = — . (34) 

n 1920 3 

We name it "EA-3" and show the performance of the three BP 

decoders when applied to EA-3 over depolarizing channels in 

Fig. [7] 




Standard BP decoder 
PC08 

Enhanced Feedback 



10'' • 

0.17 0.18 0.19 0.20 0.21 

Cross-over probability of depolarizing channel, p 

Fig. 7. Same as for Fig. 4 but using Conventional replaced by EA-3 and 
changing the maximal number of entries traversed by enhanced feedback 
decoder round down to one fortieth of the code length that is 43. 



C. Time efficiency simulation 

We evaluate and compare the time consumed by each 
decoder according to the average number of iterations 



ANoI := 



rib 



(35) 



consumed by one codeword, with n t the total number of 
iterations to finish decoding n& blocks. In our simulation 
Hi, = 20000. In Figs. 8 and 9 we show the time efficiency 
(complexity) comparison for the three decoders applied to 
Conventional, EA-1, EA-2 and EA-3. Evidently our algorithm 
has a smaller ANoI than PC08, which shows improved decod- 
ing efficiency when compared to PC08. 

D. Discussion 

From Figs. |4]|9] it is evident that our strategy yields a sig- 
nificantly lower BER with a lower ANoI for both conventional 
sparse quantum codes and EA sparse quantum codes. As an 
example of the effectiveness of our new decoder, consider the 
case p ~ 0.015, Fig. 4 shows that the enhanced-feedback BP 
decoder yields a 6 dB gain over PC08 and an 8 dB gain over 
the standard BP decoder while reducing the ANoI from 3.1 for 
PC08 to 2.2 for our enhanced-feedback BP decoder as seen 
in Fig. 8. 

For another point of view, in Fig. [5] we can see that 
to guarantee EA-1 having a decoding error below 10~ 3 , the 
standard BP decoder allows p =s 0.020 and PC08 allows p w 
0.022. This value is increased dramatically to p « 0.029 



-A— 
-O- 

-v- 



- Standard BP (Conventional) 
■ PC08 (Conventional) 

- Enhanced (Conventional) 
Standard BP (EA-1) 
PC08 (EA-1) 
Enhanced (EA-1) 




0.015 0.02 0.03 0.04 0.05 0.06 

Cross-over probability of depolarizing channel, p 

Fig. 8. Same as for Fig. 4 but replacing block error rate by average number 
of iterations and adding EA-1 as another example. 



— v- 

— A— 
-O- 

-v- 



Standard BP (EA-2) 
PC08 (EA-2) 
Enhanced (EA-2) 
Standard BP (EA-3) 
PC08 (EA-3) 
Enhanced (EA-3) 




0.09 0.11 0.13 0.15 0.17 0.19 0.21 
Cross-over probability of depolarizing channel, p 

Fig. 9. Same as for Fig. 8 but replacing Conventional and EA-1 by EA-2 
and EA-3 and changing the maximal number of entries traversed by enhanced 
feedback decoder from rounding down to one fifth of the code length to one 
tenth of the code length for EA-2 and to one fortieth of the code length for 
EA-3. 



using our enhanced feedback BP decoder with a similar ANoI 
performance as PC08. As for the decoding of EA-2 (Fig. [6]), 
when p = 0.09, the BER performance of PC08's method is 
1.75 times better than that of the standard BP decoder with 
2.78 times more ANoI, and the BER performance of our 
approach is approximately 10 times better than that of the 
standard BP decoder with only 1.3 times more ANoI. 

Notably, interesting results appear in the decoding of EA-3. 
From Fig. [7] and [9] we see that PC08 could not significantly 
improve the standard BP decoder even with high ANoI, yet our 
approach does improve the performance of EA-3 significantly 
with almost the same ANoI as the standard BP decoder 
excepting the last two nodes. 
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Our algorithm is capable of turning the detectable errors 
into new outputs whose syndromes are identical with corre- 
sponding observed syndromes. We refer to these outputs as 
valid outputs can be classified into three cases: (1) E out = E; 
(2) £ out + E but E oaX x E e S; (3) E oat x E £ S. The 
first two cases will yield correct outputs that can recover the 
sent quantum states successfully. However, it is hard to check 
whether £ out x E e S which makes distinguishing cases (2) 
and (3) difficult. Therefore, we choose (1) as the success 
criterion for the decoding result in our numerical simulation. 

Some detectable errors are turned into undetectable errors 
by our decoder. In order to show how often our algorithm re- 
sults in an undetectable error, we classify all the results yielded 
by our algorithm (not just valid outputs) when the original BP 
protocol fail^l into three cases: (1); cases (2) or (3); remain 
as detectable errors. In Fig. [10] we take Conventional as an 
example to depict the number of blocks that fall into each 
whereas the original BP protocol fails 50 times. It is evident 



50 
,_ 40 







^■Case (1) 

I I Cases (2) or (3) 

I | Remain as detectable errors 




ll 


1 




1 




j 





0.0150.02 0.03 0.04 0.05 
Cross-over probability of depolarizing channel, p 

Fig. 10. Block numbers fall into the three cases mentioned in context vs 
cross-over probability of depolarizing channel p for our enhanced BP decoder 
when applied to Conventional while the original BP decoder fails 50 times. 
The maximum number of iterations is 90. The maximal number of iterations 
between each perturbation is 40. The maximal number of entries traversed by 
enhanced feedback decoder round down to one fifth of the code length. 

to see from Fig. [I0]that most of the detectable errors will be 
turned to correct output but not undetectable errors by using 
our enhanced feedback BP decoder. 

VI. Conclusion 

We have developed an enhanced feedback BP decoding 
algorithm whose feedback adjustment strategy is not based 
solely on the syndrome but also on the channel model and on 
the individual values of the entries of the frustrated checks. 
Our approach retains the capability of breaking the symmetric 
degeneracy while also feeding back extra useful information to 
the BP decoder. Therefore, our feedback adjustment strategy 
yields a better error correcting capacity with relative less 
iterations compared to the proposed decoding algorithms for 
sparse quantum codes. We have considered three cases: the 

9 The 'fails' means that the syndrome of the output of the original BP 
decoder is not identical with the observed syndrome. In other words, a detected 
errors yielded by the original BP decoder. 



standard BP decoder adapted from classical decoding to de- 
code quantum codes, the superior BP decoder with Poulin and 
Chung's random perturbation IflOl , and the BP decoder with 
our new feedback adjustment strategy introduced here. 

We used the GF{A) representation of quantum codes to 
construct our feedback-based algorithm and exploit not just 
the syndrome but all the measurement resulting from stabilizer 
measurements. Then we used the block error rate (BER) 
and average number of iterations (ANoI) to demonstrate a 
dramatic better error correcting capacity improvement with 
relative less iterations. This result was achieved using new 
feedback adjustment strategy vs the two alternatives: standard 
BP decoder and the BP decoder with random perturbation. 

Appendix 

As shown in Section V, the net rate for EA-2 and EA-3 
are and -1/3 respectively. As pointed out by Brun, Devetak 
and Hsieh, "in general, net rates for EA quantum codes can be 
positive, negative, or zero. The zero rate does not mean that 
no qubits are transmitted by this code! Rather, it implies that 
a number of bits of entanglement is needed that is equal to 
the number of bits transmitted." Compared to quantum codes 
with high net rates, EA quantum codes with zero or negative 
net rates employ many more physical qubits including a great 
number of bits of entanglement that being assumed to be error- 
free0 thus these codes can tolerate stronger noise. Therefore, 
it is not surprising that EA-2 and EA-3 greatly increase p the 
cross-over probability of depolarizing channel when compared 
to quantum codes with similar code length and higher net rates. 

For comparison, based on "Construction B" |9j (same as the 
construction of Conventional), we construct a quantum code 
with a similar code length as EA-2 but having a relative higher 
net rate. We name this example as EX which has a code length 
of 800 and a net rate 1/2. In Fig. QT| we see that EX has 
comparable BER performance to other examples presented in 
existing literatures which have similar parameters as EX but 
significantly worse than EA-20 
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